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Abstract 


Let  P  be  a  simple  polygon  with  N  vertices,  each  assigned  a  weight  c;  (O^i $N)  with  c.€{0,l}.  We  define  die 
weight  C  of  P  as  the  added  weight  of  its  vertices,  Le.,  C  =  c^+...+cff.  Making  the  assumption  that  die 
vertices  of  P  have  been  sorted  along  some  axis  --which  can  be  done  in  0(Nlog  N)  time--,  we  prove  that  it  is 
possible,  in  O (JV)  time,  to  find  two  vertices  a,b  in  P,  such  that  the  segment  ab  lies  entirely  inside  the  polygon  P 
and  partitions  it  into  two  polygons,  each  with  a  weight  not  exceeding  2C/3.  We  also  give  a  list  of  problems 
which  can  be  solved  efficiently  with  a  divide-and-conquer  strategy  based  on  that  result 
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1.  Introduction 

It  is  a  fact  that  fast  algorithms  often  owe  their  efficiency  to  a  handful  of  mathematical  truths  that  bring  to 
light  unsuspected  specificity  ~  or  singularity  --  about  the  problem  under  consideration.  But  besides 
recognizing  their  mere  aptitude  at  serving  computational  purposes,  one  may  wish  to  add  a  bit  of  classification 
into  the  nature  of  those  subservient  truths. 

To  ease  the  matter,  let  us  take  two  examples:  proving  that  a  minimum  convex  decomposition  of  a  polygon 
could  always  be  found  free  of  interior  parts  was  a  major  step  in  deriving  a  polynomial  algorithm  for  that 
problem  [CH80].  Similarly,  establishing  the  unimodality  of  the  vertex-coordinate  function  of  a  convex 
polygon  was  instrumental  in  the  setting  of  logarithmic  intersection-algorithms  for  convex  objects  [CD80]. 
There  is,  however,  a  notable  difference  in  methodology  between  these  two  cases.  Whereas  the  first  example 
presents  us  with  a  mathematical  fact  which  seems  inherently  helpful  for  deriving  any  algorithm,  die  latter 
brings  out  a  compelling  flavor  of  binary  search-like  technique  which  strongly  suggests  the  proper  algorithmic 
treatment 

Another  case  which  witnesses  both  phenomena  is  the  well-known  nearneighbor  problem,  in  which  all 
nearest  neighbor  pairs  in  a  given  set  of  points  are  to  be  computed.  One  of  the  most  efficient  algorithms  for 
this  problem  relies  on  an  elegant  geometric  construction  known  as  the  Voronoi  diagram  [SH77].  Although  the 
existence  of  this  diagram  was  undoubtedly  die  keystone  in  the  elaboration  of  the  algorithm,  it  fell  short  of 
even  suggesting  an  effective  method,  since  there  was  no  immediate  evidence  that  constructing  the  Voronoi 
diagram  was  to  be  any  easier  than  computing  the  neighbor  pairs  directly.  Another  mathematical  fact,  Le.,  the 
decomposability  of  the  Voronoi  construction,  was  indeed  needed  to  fire  the  final  blow  and  crack  the  problem. 
There  we  should  observe  that  although  both  facts  are  geometric  in  nature,  the  latter  is  combinatorial  in  spirit, 
as  it  draws  its  motivation  from  the  algorithmic,  all-purpose,  divide-and-conquer  technique. 

Lipton  and  Taijan’s  planar  separator  theorem  [LT77.L177]  is  a  notable  example  of  a  systematic  technique 
for  introducing  a  computational  tool,  Le.,  divide-and-conquer,  into  a  whole  class  of  related 'problems,  Le., 
planar  graph  problems.  Drawing  its  inspiration  from  this  philosophy,  this  paper  presents  a  theoretical  result  on 
polygon  decomposition  which  can  be  applied  to  derive  a  number  of  efficient  algorithms  for  geometric 
problems,  in  particular,  problems  of  convex  decompositions,  triangulation,  visibility,  and  internal  distance. 

2.  The  Polygon-cutting  Theorem 

Let  P  be  a  simple1  polygon  with  vertices  v1,...,vN  in  clockwise  order.  Let  OXY  be  an  orthogonal  system  of 
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A  polygon  is  said  to  be  simple  iff  only  adjacent  edges  intersect 
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reference.  Wlog,  we  can  always  assume  that  no  two  vertices  have  the  same  X-coordinate.  We  define  a  partial, 
so-called  vertical,  order  as  follows: 

We  say  that  "edge  e  >  edge  f  iff  their  projections  onto  the  .Y-axis  overlap  and,  restricting 
ourselves  to  the  overlapping  area,  e  lies  entirely  above  /(fig.l). 

[  FIGURE  1 1 

Figure  1:  The  partial  order  among  the  edges  of  P. 

Throughout  this  paper,  we  will  assume  that  along  with  a  description  of  the  boundary  of  P,  provided  by  a 
doubly-linked  list  LP,  we  have  available  both  a  doubly-linked  list  LV  of  all  the  edges  in  topological  ( vertical, ) 
order,  and  a  doubly-linked  list  LH  of  all  the  vertices  sorted  by  Y-coordinates  [KN73J.  The  preprocessing 
involved  in  setting  up  these  lists  requires  0(Mog  N)  time  and  O (N)  space.  The  decomposition  algorithm 
which  we  will  describe  later  on  runs  in  linear  time,  with  this  preprocessing  in  hand.  Note  that  it  is  legitimate  to 
separate  both  tasks,  since  in  the  applications  which  we  will  mention,  the  decomposition  algorithm  will  be 
called  recursively  several  times,  while  the  preprocessing  will  be  needed  initially,  once  and  for  all 

The  goal  of  this  paper  is  to  prove  the  following  theorem: 


Theorem  1:  The  Polygon-cutting  Theorem.  Let  P  be  a  simple  polygon  with  N  vertices  v^...,^ 
each  assigned  a  weight  c{  (c.=0,l).  Let  C(P)  denote  the  total  weight  of  P,  defined  as  the  sum 
Cj + ... + Cjp  and  assume  that  C(P)>Z  With  die  lists  LP.LH.LV in  hand,  it  is  possible  to  find,  in 
0(N)  time,  a  pair  of  vertices  Vj,v.  such  that  the  segment  vv.  lies  entirely  inside  the  polygon  P  and 
partitions  it  into  two  simple  polygons  PVP2  satisfying: 

CiPJsdPJzlCim 

The  weights  of  the  vertices  in  Px  and  P2  are  die  same  as  in  P,  except  for  v.  and  v.,  for  which  we  will  assume 
that  in  both  Px  and  P2,  these  weights  become  0.  This  assumption  is  made  only  for  the  sake  of  simplicity,  and 
other  conventions  (e.g„  keeping  the  same  weights  Cj.c j  in  both  P2  and  /*2)  are  indeed  acceptable,  if  we  are 
ready  to  add  a  term  +2  to  2C(f)/3  in  the  inequality  of  Theorem  1.  To  facilitate  our  task,  we  will  first  prove 
the  theorem  with  slighdy  relaxed  requirements. 


2.1 .  An  existence  theorem 

To  begin  with,  we  will  prove  the  existence,  not  of  two  vertices,  but  of  two  points  on  the  boundary  of  P, 
satisfying  the  inequalities  of  Theorem  1. 

Theorem  2:  Same  assumptions  as  Theorem  1.  There  exists  a  pair  of  points  A.B  on  the 
boundary  of  P,  such  that  the  segment  AB  is  parallel  to  the  Y-axis,  lies  entirely  in  the  polygon  P, 
and  partitions  it  into  two  simple  polygons  PyP2  satisfying: 

C(P2)  £  C(P2)  5  2C(P)/3 


If  A  and  5  are  not  vertices  of  P,  for  consistency,  we  assign  them  a  0-weight  The  underlying  notion  of 
"equal  size”  decomposition  expressed  in  the  polygon-cutting  theorem  motivates  the  introduction  of  a  distance 
fiinction  d (A,B),  defined  between  two  points  A,B  on  the  boundary  of  P  as  the  minimum  path  weight  between 
A  and  B.  More  precisely,  let  v.v.+1  (resp.  v.v.+1)  be  the  edge  of  P  containing  A  (resp.  B).  If  A  (resp.  B)  is  a 
vertex  of  P,  it  is  assumed  to  be  v.  (resp.  vp.  We  introduce  the  function  h,  defined  as  follows2: 

h(.d,5)  =  <?j+i+tj+2+  -  +ci 

from  which  we  can  define  d 1(44): 

dUfl)  =  minlh(y4,23),C(^)  -  hM,*)] 

Starting  at  the  edge  v^,  we  label  each  edge  of  P  recursively,  as  follows: 

XfVjVj)  =  cx 

X(vivi+i)  =  X(VHvp  -I-  cs 

Note  that  this  labeling  gives  us  an  alternate  way  of  defining  the  distance  between  two  boundary  points  A,B: 
d (A,B)  m  mini |X(VjV.+1)  -  MVjVj+1)| , C(/')-|X(vivi+1)  -  Myj+1)| I 

We  are  now  in  a  position  first  to  prove  the  existence  of  the  segment  AB,  as  defined  in  Theorem  2,  then  to 
describe  an  efficient  method  for  finding  it  As  we  will  see,  the  first  step  is  not  superfluous;  it  is  an  essential 
ingredient  in  ensuring  the  correctness  of  the  algorithm. 

Choose  the  leftmost  point  of  P  as  the  starting  point  of  the  left-to-right  sweep  of  a  vertical  segment  S=AB 
{A  below  B).  S  will  always  stretch  vertically  so  as  to  keep  its  endpoints  A,B  in  permanent  contact  with  foe 
boundary  of  P.  It  will  thus  be  able  to  move  continuously  to  the  right,  until  it  must  either  expand  (fig.2.1)  or 
split  (fig.2.2).  At  any  time  during  the  course  of  the  motion,  S  will  be  assigned  a  value  A=h(/4,5)  to  indicate 
how  close  it  is  to  being  the  desired  segment  We  observe  that  initially,  A=C (P),  and  that  as  long  as  S  moves 
continuously,  A  decreases  monotonously  by  unit  steps.  When  either  situation  depicted  in  Figure  2  arises,  we 
can  always  write 

A  =  h  (A,B)  =  hUQ  +  h(C,5)  (1) 

from  which  we  can  derive  a  decision  procedure  for  redefining  S. 

Starting  from  the  leftmost  vertex  of  /*,  move  S  from  left  to  right  stretching  or  shrinking  this 
segment  so  that  it  entirely  lies  in  P,  and  its  endpoints  always  lie  on  the  boundary  of  P.  As  long  as 
the  motion  of  S  is  continuous,  check  whether  A>2C(/>)/3.  in  which  case  continue,  else  stop. 
When  falling  in  cither  case  of  fig.2,  reset  S  to  AC  if  h(/f,C)>h(C,/?),  or  to  CB  otherwise.  Note  that 
if  S  is  reset  to  CB  in  the  case  of  fig.2.1,  the  motion  must  reverse  its  direction. 
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Relation  (1)  shows  that  every  discontinuity  causes  A  to  decrease  by  at  most  half,  while  otherwise  A 
decreases  at  most  by  unit  steps,  since  we  have  assumed  that  no  two  vertices  may  lie  on  the  same  vertical  line. 
Since  A  eventually  vanishes  and  C(P)>2,  it  must  take  on  some  value  in  the  interval  [C(f)/3,2C(/1)/3],  at 
which  point  the  procedure  will  stop  and  return  the  desired  segment  AB.  This  completes  the  proof  of 
Theorem  1  □ 

{ FIGURE  2  J 

Figure  2:  Proving  the  existence  of  a  separator  S. 


2.2.  A  relaxed  version  of  the  polygon-cutting  theorem 

Unfortunately,  Theorem  2  falls  short  of  providing  an  efficient  algorithm  for  computing  AB.  We  can, 
however,  graft  to  it  a  binary  search-like  structure  to  improve  the  performance  of  a  naive  implementation.  The 
purpose  of  this  section  is  thus  to  prove  the  following  result: 

Theorem  3:  Same  assumptions  as  Theorem  1.  It  is  possible,  in  O (N)  time,  to  find  a  pair  of 
points  A.B  on  the  boundary  of  P,  such  that  the  segment  AB  lies  entirely  in  the  polygon  P,  and 
partitions  it  into  two  simple  polygons  PVP2  satisfying: 

CiPj)  S  CfPp  s  2C(/)/3 

The  algorithm  which  we  will  describe  in  order  to  prove  Theorem  3  is  recursive;  it  requires  0(N)  time  to  cut 
down  the  size  of  the  problem  by  half,  therefore  its  overall  performance  is  linear. 

1.  Recall  that  initially,  we  have  available  both  the  vertical  topological  order  (LV)  of  die  edges  of  P 
and  the  horizontal  order  (LH)  of  its  vertices.  The  latter  list  permits  us  to  determine  a  median 
vertical  line  L  in  O(A0  time,  i.e„  a  line  parallel  to  the  K-axis  that  separates  the  vertices  of  P  into 
two  sets  of  size  \CXP)/1\  and  [C(/J)/2j  respectively.  Once  again,  this  is  always  possible  since  no 
two  vertices  may  have  the  same  X-coordinate  (fig.3). 

2:  Compute  all  the  intersection  points  of  L  with  the  boundary  of  P,  and  sort  them  by  F-coordinates, 
doing  all  of  this  in  O (N)  time  with  the  list  LV.  Next,  form  the  intersection  segments  AxBv„.,AtBx 
in  ascending  order  (fig.3). 

3.  If  M  is  a  point  lying  on  an  edge  VjVj+]  of  P,  distinct  from  its  endpoints,  by  extension,  we  define 
X(A/)  as  X(VjVi+ 1).  Check  whether  any  of  the  segments  A.B.^  satisfies  the  relation 

C(P)/3  <  |X(/1.)  -  X(0;)|  <  2C(F)/3,  for  i=l,...Jt. 

If  yes,  the  corresponding  segment  A.XB.X  can  be  chosen  as  AB ,  and  the  procedure  can  stop. 
Otherwise,  go  to  4. 

4.  Transform  the  list  LP  as  follows.  For  every  A.xB~x  (i= l,..,k)  in  turn,  if  |X(/f.)  -  X(0j)|  <  C(/V3, 
start  at  A.  or  B.  (say,  AJ  whichever  has  the  smaller  label  X(/l)  or  X(^),  and  proceed  to  traverse 
LP  in  clockwise  order,  until  B.  is  reached,  deleting  all  the  vertices  and  links  visited.  Finally,  close 
the  chain  by  inserting  A.  and  0.  into  LP,  setting  a  double  link  between  the  two  points.  If 
|X(/fj)  -  X(0.)|  >  2C(/>)/3,  perform  the  same  sequence  of  operations,  now  starting  at  the  point  Ai 
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or  B{  with  the  larger  label  X(/4j)  or  X(R). 

5.  Using  the  new  description  LP  of  the  pruned  polygon,  update  the  lists  L  V  and  LH  by  deleting  the  I 

edges  and  vertices  no  longer  existing  in  LP. 

6.  Apply  the  algorithm  recursively  until  termination. 

[FIGURE  31 

Figure  3:  The  recursive  algorithm  for  computing  the  separator  S.  * 

A  few  words  of  explanation  may  be  needed  to  justify  the  algorithm.  To  be  begin  with,  observe  that  once  die 
intersection  points  are  available  in  vertical  order,  setting  up  the  list  of  segments  AiBi  is  straightforward,  since 
the  list  of  consecutive  segments  formed  by  the  intersection  points  alternates  inside  and  outside  P.  Also,  it  is  ! 

easy  to  verify  that  the  deletion  of  edges  is  accompanied  by  the  setting  of  links  which,  in  particular,  ensures 
that  no  edge  is  visited  more  than  twice.  It  may  be  the  case  that  handling  a  segment  A^  may  also  handle 
several  others  at  the  same  time.  For  example,  see  the  effect  in  fig.4  of  handling  AjB^  before  AaB4,  or 
vice-versa.  In  this  figure,  the  dashed  area  is  to  be  removed  from  P.  Finally,  we  note  that  deletions  are  indeed  I 

permissible  since  any  part  of  P  cut  off  by  A{B{  weighs  at  most  d (A^B^  <  C (P)/3.  Thus  any  vertical  segment 
AB  in  that  part  has  a  distance  d(A,B)  <  C (P)/3,  and  may  therefore  be  dismissed.  Since,  when  iterating  down 
the  recursion  level,  all  the  segments  Afi  are  being  cut  off  and  their  attached  part  removed,  there  remains  only 
one  connected  part  which  lies  entirely  on  one  side  of  the  median  line,  and  to  which,  therefore,  P  can  be  reset  I 

i 

before  calling  the  algorithm  recursively.  Of  course,  the  labeling  X  remains  unchanged  The  correctness  of  die 
algorithm  follows  from  Lemma  2  and  from  the  fact  that  we  remove  only  non-candidates.  The  reason  for 
choosing  a  median  line  to  operate  the  cut-offs  is  now  apparent  Since  all  the  vertices  on  one  side  of  the  line  are 
bound  to  disappear,  the  recursion  will  be  invoked  on  a  problem  of  size  reduced  by  at  least  half.  This  I 

completes  the  proof  of  Theorem  3.  □ 

[FIGURE  4J 

Figure  4:  Pruning  the  polygon  P. 

2.3.  Completing  the  proof  of  the  polygon-cutting  theorem 

We  may  now  turn  our  attention  back  to  Theorem  1.  Let  VjVj+1  and  y.+1  be  the  edges  of  P  that  contain 
the  points  A  and  B  of  Theorem  3,  respectively.  To  prove  the  desired  result,  one  may  be  tempted  to  slide  A  and 
B  towards  the  endpoints  of  v.v.+1  and  y.  +  1  respectively,  until  one  of  the  configurations  y,  v.Vj+1,  vi+JVj,  or 
v. + 1  v. + x  has  been  reached.  Unfortunately,  obstacles  may  prevent  this  from  ever  happening  (fig.5),  so  our  next 
step  will  be  to  take  a  closer  look  at  these  possible  obstacles. 

Since  the  quadrilateral  yivi+ivj vj+lvi  conta*ns  ^e  segment  AB,  it  is  a  simple  polygon,  and  AB  partitions  it 
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into  two  polygons  Q1  =  Av.^^BA  and  Q1  =  ABv.  ^A.  As  a  corollary  of  the  Jordan  Curve  Theorem 
[HSSS]3,  it  is  easy  to  show  that  the  only  obstacles  encountered  in  Ql  (resp.  Q2)  are  vertices  in  the  set 

{vi+i*vi+2’-’vj}  (resP-  {vj+rvj+2 . vi^‘  Moreover,  the  segment  S=AB  can  encounter  only  vertices  on  the 

convex  hull  Hx  (resp.  H2)  of  the  vertices  of  P  lying  inside  Qx  (resp.  Q2)  (fig.5),  as  is  shown  in  following  result 

Lemma  4:  The  segment  AB  intersects  any  edge  of  P  (outside  of  A  or  B)  if  and  only  if  it 
intersects  the  boundary  of  either  Hx  or  Hr 

Proof:  Since  Hx  (resp.  //,)  lies  entirely  inside  Qx  (resp.  QX  AB  intersects  any  edge  of  P  lying 
in  Q  iff  the  infinite  line  passing  through  AB  does,  hence  iff  AB  lies  outside  of  Hx  and  Hr  □ 

The  next  task  is  to  compute  the  convex  hulls  and  Hr  We  only  give  the  details  of  the  algorithm  for  Hy 
the  other  case  being  strictly  similar.  Since  we  cannot  afford  to  use  a  standard  0(Mog  N)  algorithm  to  simply 
compute  the  convex  hull  of  the  vertices  of  P  in  Qy  we  must  exploit  the  fact  that  these  vertices  lie  on  a 
polygonal  line  in  order  to  achieve  linear  time. 

[FIGURE  51 

Figure  5:  Defining  the  domain  of  safety  for  AB. 

To  begin  with,  let  us  give  an  informal  description  of  the  algorithm.  The  goal  is,  in  a  first  stage,  to  produce  a 
polygon  K2  which  lies  entirely  in  Qv  and  whose  convex  hull  is  exactly  Hy  Kx  consists  essentially  of  polygonal 
chains4  L  made  of  consecutive  edges  from  the  set: 

L  =  {  vj+1vi+2 ,  vj+2vi+3  ......  v>1vj } 

Each  chain  has  the  property  that  it  lies  in  Qx  and  intersects  v.+1v.  in  two  points,  Uk  and  Kk.  Moreover,  no 
two  segments  t/kKk  overlap  (fig.6.1).  To  compute  these  chains,  we  must  distinguish  between  two  types  of 
edges  in  L.  An  edge  vkvk+1  is  said  to  be  entering  (resp.  exiting)  if  it  intersects  vj+1Vj,  and  vk+1  lies  inside 
(resp.  outside)  Qy  The  algorithm  proceeds  as  follows: 

Wlog,  assume  that  there  is  at  least  one  vertex  vj+2  from  vj+1  to  v.  in  clockwise  order.  Traverse  L  from 
vi+lvi+2 t0  vj-iyj’  st0PP'n8  at  entering  and  exiting  edges  and  taking  the  following  actions.  If  the  current  edge 
is  entering,  it  may  be  the  endpoint  £/k  of  a  new  chain  L^.  To  decide  of  this,  look  at  the  next  exiting  edge  in  L; 
if  it  intersects  vj+1v.  in  a  point  Vk  on  the  segment  C/kVj,  we  have  indeed  a  new  chain  Lk  from  Uk  to  Kk. 
Otherwise,  not  only  don’t  we  have  a  new  chain  Lk,  but  the  chain  just  visited  may  enclose  previously  computed 
chains,  which  must  then  be  deleted.  For  that  purpose,  we  use  a  stack  to  hold  the  pairs  (f/kKk  ,  Z^),  so  that 

^Th  e  Jordan  Curve  Theorem  states  that  a  dosed  curve  in  the  plane  partitions  the  plane  into  two  connected  regions:  the  inside  and  the 
outside. 

A 

We  define  a  polygonal  chain  as  any  connected  sequence  of  segments  such  that  only  consecutive  segments  may  intersect 
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deletions  may  be  done  efficiently.  The  algorithm  is  straightforward,  so  we  may  omit  the  details. 

Initially,  the  stack  is  empty,  and  the  current  edge  e  is  the  first 
entering  edge  of  L 

begin 

Let  /be  the  next  exiting  edge  in  L  following  e,  and  let  U,V  be 
respectively  the  intersections  of  e  and  /with  vi+1Vj. 

if  Flies  on  Uv. 
then 

Let  Lu  be  the  chain  between  e  and  fin  L. 

Push  (UV,  LJ  onto  stack. 

Go  to  next  entering  edge  e'  in  L,  whose  intersection 
with  vi+1v.  lies  on  Kv.+1,  then  iterate. 

else 

Go  to  next  entering  edge  e  whose  intersection  If 
with  v.+ jV.  lies  on  vj+ J V. 

Pop  aU  pairs  (Ut  Kk)  off  die  stack  as  long  as  Kk  lies  on  U"Vy 
Iterate, 
end 

[  FIGURE  6  ] 

Figure  6:  Computing  the  polygon  Kt 

To  prove  the  correctness  of  the  algorithm,  we  begin  by  observing  that  the  intersection  of  L  with  Qx  consists 
of  chains  whose  endpoints  lie  on  v.+  x  Vj,  and  that  Kx  consists  of  exaedy  all  the  maximal  chains.  A  chain  is  said 
to  be  maximal  if  it  does  not  lie  in  the  enclosure  of  any  other  chain  with  v.+1Vj.  A  maximal  chain  is  also 
characterized  by  the  fact  that  the  segment  formed  by  its  endpoints  does  not  lie  inside  any  other  such 
segments. 

From  the  Jordan  Curve  Theorem  [HS55],  we  derive  the  fact  that  a  maximal  chain  from  U  to  V,  in  clockwise 
order,  has  its  endpoint  Flying  above  U (i.e.,  on  the  segment  £Ap.  In  consequence,  only  the  chains  which  move 
towards  *  are  candidates  for  being  part  of  Kv  which  justifies  the  selection  criterion  of  the  algorithm.  On  the 
other  hand,  we  can  also  show  that  a  non-maxima]  chain  which  moves  towards  v.  is  necessarily  enclosed  by 
another  chain  moving  away  from  Vj.  This  explains  the  deletion  rule.  Finally,  the  last  observation  to  make  is 
that  a  chain  from  U  to  V  which  moves  away  from  v.  (i.e.,  U  lies  on  Kv.)  must  be  enclosed  by  a  subsequent 
maximal  chain,  therefore  since  maximal  chains  are  computed  "towards”  v,  we  may  skip  directly  to  the  next 
entering  edge  in  L  that  intersects  v.+1v.  at  a  point  IT  below  K(i.c.,  If  lies  on  v.+1K)  -  sec  illustration  of  the 
various  cases  in  fig.6.1.  This  completes  the  proof  of  correctness,  and  shows  that  all  the  chains  may  be 


computed  in  sorted  order  along  the  segment  vj+ ^  all  these  computations  requiring  0(N)  time. 

The  final  step  in  computing  Kx  is  to  connect  all  the  chains  Lk  together  in  the  order  in  which  they  appear  in 
the  stack.  To  do  so,  we  borrow  segments  from  vi+1v.,  as  shown  in  fig.6.2.  We  may  now  apply  any  standard 
linear  convex  hull  algorithm  for  simple  polygons  [LE80.SH77]  in  order  to  obtain  Hx  in  0(A)  time. 

Assuming  that  both  H1  and  H2  are  available,  we  are  now  in  a  position  to  give  an  algorithm  for  finding  the 
two  vertices  of  Theorem  1.  The  idea  is  to  connect  the  vertices  of  H1  with  those  of  H2,  so  as  to  triangulate  the 
polygon  tf  defined  as  the  area  between  Hl  and  H2  containing  AB.  We  claim  that  at  least  one  of  the  edges  of 
the  triangulation  will  provide  the  desired  pair  of  vertices,  with  the  property  of  Theorem  1.  The  algorithm 
proceeds  as  follows: 

Let  /^....,/ip  and  be  the  vertices  of  Hl  and  H2,  respectively,  as  we  traverse  them  from  v.+1vt  to 

VjVj+t*  ‘-e-’  ^i= vi+i*  hp~  vj*  *i=  vi*  *q=  vj+ r  We  maintain  two  pointers,  h  on  Hl  and  k  on  H2,  moving  them 
from  Aj  to  Ap  and  to  Jfcq,  respectively,  and  computing  the  triangulation  on  the  fly.  Note  that,  at  all  times,  the 
segment  hk  intersects  H,  and  H2  only  at  its  endpoints  (fig.7.1). 

[FIGURE  7] 

Figure  7:  Triangulating  the  domain  of  safety. 

The  simplest  way  of  describing  the  algorithm  is  recursively.  Initially,  hk  is  v^v;  the  algorithm  terminates 
with  hk  =  VjVj+1. 

Let  hk  =  At*u,  and  consider  the  quadrilateral5  hht+  xku+1kh.  We  will  show  in  Lemma  5  that  at 
least  one  of  its  diagonals,  hku +1  or  kh{ + j,  connects  Hl  and  H2  without  intersecting  these  polygons 
outside  of  its  endpoints,  i.e.,  lies  entirely  in  H  .  Moreover,  this  diagonal  can  be  found  in  constant 
time.  We  may  then  determine  that  diagonal,  add  it  to  the  triangulaiion,  set  hk  to  it,  and  iterate. 

The  algorithm  clearly  runs  in  linear  time.  Also,  the  assurance  that  it  effectively  produces  a  triangulation  of 
H*  comes  from  the  fact  that  it  keeps  only  edges  which  lie  entirely  in  //*,  and  that  the  pointers  h  and  k  pass  a 
vertex  only  after  a  diagonal  has  been  assigned  to  it  Thus  there  only  remains  to  prove  the  following  lemma: 

Lemma  5:  If  the  segment  h^  connects  Hx  and  H2  and  lies  entirely  inside  H *,  so  does  one  of 
the  diagonals  \k^Jrl  or  kji^  v  Moreover,  this  diagonal  can  be  found  in  constant  time. 

Proof:  Consider  the  line  passing  through  hk,  oriented  from  h  to  k.  If  a  point  lies  to  the  right 
(resp.  left)  of  this  line,  we  will  say  that  it  lies  below  ( resp.  above )  hk.  Since  Hx  and  H2  are  convex,  at 

5For  consistency,  we  define  +  j  and  +  ^  as  v.  +  ^  and  Vj,  respectively. 
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least  one  of  the  vertices  vj  or  v.+ 1  lies  above  hk.  therefore  it  is  impossible  that  both  h{+l  and  &u+1 
lie  oelow  hk.  Indeed  this  would  involve  the  existence  of  at  least  three  intersection  points  between  a 
line  and  a  convex  boundary,  leading  to  a  contradiction.  If  only  one  segment,  say  A[+1  lies  above 
hk,  it  can  be  easily  determined  in  constant  time,  and  since  in  that  case,  all  of  H2  lies  below  hk,  the 
diagonal  \+lk  does  not  intersect  H2  (nor  Hl  either)  outside  of  its  endpoints,  and  may  thus  be 
chosen  as  the  next  segment  of  the  triangulation  (fig.7.2).  If,  on  the  other  hand,  both  h  2  and  *u+, 
lie  above  hk,  the  quadrilateral  hh{+  xku  +  x  k  is  a  simple  polygon,  therefore  it  contains  at  least  one  of 
its  diagonals  entirely  (fig.7.3),  and  this  diagonal  can  be  found  in  constant  time.  Note  that,  because 
of  its  convexity,  Hl  (resp.  H2)  lies  totally  on  one  side  of  the  line  passing  through  hht+l 
(resp.  **u+1),  therefore  the  whole  quadrilateral,  hence  the  chosen  diagonal,  lies  inside  the  polygon 
H*,  which  completes  the  proof.  □ 


The  purpose  of  triangulating  the  polygon  H*  will  become  apparent  with  the  following  result 


Lemma  6:  There  exists  an  edge  uv  in  the  triangulation  of  H*  which  satisfies  the 
relation:  C(P)/3  <  h(u,v)  $  2C(P)/3. 


Proof:  From  Theorem  3,  we  know  that  AB  partitions  P  into  two  polygons  with  weights  between 
C(P)/ 3  and  2C(P)/3,  which  gives  the  relations 

C(P)/3  <  min  ( h(A,B),h(B,A) )  <  max(  h(A,B),h(B,A)  )  <,  2C(P)/3 

As  a  result  any  pair  of  vertices  cub  on  Hl  (resp.  H^),  with  b  following  (resp.  preceding)  a  in  the 
list  {hv...,h  }  (resp.  {fcj, ...,£})  satisfies  the  relation: 

h (b,a)  <  2C(P)/3  (1) 

On  the  other  hand,  each  triangle  abc  of  the  triangulation  has  one  side  ab  on  the  boundary  of 
either  Hx  or  H2,  with  the  two  others  ac,  be  constructed  by  the  triangulation  algorithm.  ''Hog,  let  ac 
be  the  segment  of  the  triangle  constructed  first  (i.e.,  ac  lies  below  be).  We  always  have 

h(c,a)  =  h(b,a)  +  h(c,b)  (2) 

Now  we  can  show  that  a  simple  upward  scan  through  the  faces  of  the  triangulation,  i.e.,  starting 
at  the  triangle  adjacent  to  v.v.+1  and  ending  at  the  triangle  adjacent  to  v.k+1,  will  inevitably  lead 
to  the  desired  segment  of  Theorem  1.  To  see  that,  we  may  obviously  assume  that  none  of  the  edges 
ab  of  Hx  or  H2  satisfies  the  relations: 

C(P)/3  <  h(M  <  2C(P)/3 , 
otherwise,  we  have  achieved  our  goal. 

In  that  case,  Relation  (1)  shows  that  for  any  triangle  abc  visited,  the  edge  on  the  boundary  of 
W*,  say  ab.  satisfies  the  stronger  inequality 

h(M  <  C(P)/3 , 

which,  combined  with  Relation  (2),  leads  to 

h(c,b)  >  h(c,a)  -  C(P)/3. 

Since  h(vi,vj+|)=C(/3)-cj+1  and  h(v+1,v.)=Cj+1,  it  follows  that  if  a^byC^b^  ...  is  the  sequence 
of  interior  edges  visited  in  the  traversal  of  the  triangulation,  with  the  points  am  (resp.  6m)  on  Hx 
(resp.  H2),  the  sequence  hfaj.i^.hfaj,^),...  is  monotonously  decreasing  from  C(P)— cj+1  to  Cj+, 
by  jumps  of  at  most  C(P)/3.  In  consequence,  it  must  take  on  at  least  one  value  in  the  interval 
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[C(P)/3 , 2C(P)/3],  which  can  be  chosen  as  the  pair  u,v.  □ 

The  proof  of  Theorem  1  is  now  complete.  Computing  H1  and  H2  definitely  constitutes  the  most  difficult 
part  of  the  algorithm  to  implement  We  may  observe,  however,  that  this  overhead  will  often  be  unnecessary 
since,  in  practice,  it  may  be  seldom  the  case  that  the  segment  AB  of  Theorem  3  is  prevented  from  sliding 
towards  the  endpoints  of  its  supporting  edges. 


3.  Applications  to  polygon  decomposition  problems 

It  is  intuitive  that  the  polygon-cutting  theorem  should  lead  to  efficient  methods  for  partitioning  a  polygon 
into  convex  pieces.  We  will  examine  two  instances  of  this  problem:  in  one,  what  is  desired  is  a  partition  of  the 
polygon  into  a  small  number  of  convex  pieces,  while  in  the  other,  only  a  triangulation  of  the  polygon  is 
sought,  without  consideration  of  optimality6. 

3.1 .  Convex  decompositions 

Given  a  simple,  non-convex  polygon  P,  find  a  minimum  number  of  convex,  pairwise  disjoint 
polygons,  whose  union  is  P. 

This  problem  has  been  well-studied  [CH80,CD79,FS81,GJ78,SC78,SV80,T080],  and  several  algorithms 
have  been  discovered  for  producing  minimal  or  near-minimal  decompositions.  Here  we  consider  only 
decompositions  which  do  not  introduce  new  points,  i.e.,  all  the  vertices  of  the  polygons  are  vertices  of  P. 

In  connection  with  the  previous  section,  we  will  assign  to  each  vertex  v.  of  P  a  weight  c  =  l  if  its  adjacent 
edges  form  a  reflex  angle  (in  which  case,  v.  is  called  a  notch),  and  a  weight  c4=0  otherwise.  Thus  we  can  apply 
the  polygon-cutting  theorem  (Theorem  1)  iteratively  to  decompose  P  into  smaller  polygons.  Note  that  since, 
with  our  convention,  the  endpoints  of  the  splitting  segment  lose  their  weights,  the  number  of  notches  CVC2  of 
the  two  parts  is  each  bounded  by  2C/3+2,  where  C  is  the  number  of  notches  in  the  original  polygon.  As  a 
result,  we  must  stop  the  iteration  when  the  algorithm  ceases  to  reduce  the  number  of  notches,  i.e.,  when  all  the 
parts  have  a  number  of  notches  satisfying:  C<2C/3+2,  i.e.,  C< 6.  Finally,  to  resolve  the  remaining  reflex 
angles,  we  consider  each  of  them  in  turn,  proceeding  as  follows: 

Let  Q  be  the  polygon  (with  at  most  6  reflex  angles),  and  v  be  the  notch  exhibiting  the  reflex 
angle  to  be  resolved.  Let  L  (resp.  R)  denote  the  ray  (i.e.,  semi-infinite  line)  starting  at  v  in  the 
direction  of  the  edge  ending  (resp.  starting)  at  v  (fig.8).  Compute  the  intcrscction(s)  of  L  (resp.  R) 


Applying  a  quality  criterion  to  the  triangulation  of  a  polygon  is  common  practice  in  numerical  analysis,  where  an  area  distribution  or 
a  shape  function  is  often  to  be  optimized.  There  are  many  good  reasons,  however,  for  making  the  availability  of  any  triangulation 
desirable.  Sec  [Cl  182],  for  example. 
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with  the  boundary  of  Q,  and  keep  only  the  intersection  point  A  (resp.  B )  closest  to  v.  If  A  and  B  lie 
on  different  edges,  there  exists  at  least  one  vertex  on  the  part  of  the  boundary  of  P  between  A  and 
B  which  can  be  joined  to  v,  so  as  to  resolve  the  reflex  anglft  at  v  (fig.8.1).  For  example,  we  can 
choose  the  vertex  w  between  A  and  B  that  minimizes  the  angre  (vB,vw),  while  keeping  it  positive. 

If.  on  the  other  hand,  A  and  B  lie  on  the  same  edge  vAvB  (fig.8.2),  we  compute  the  vertex  a  of  the 
list  (v,...,v^>,  given  in  clockwise  order,  which  lies  in  the  triangle  vAva  and  minimizes  the  angle 
(v^,va).  Similarly,  we  compute  the  vertex  b  which  lies  in  vBvg  and  minimizes  the  angle  (vb,vB). 

Both  of  these  operations  can  be  executed  in  linear  time.  Note  that  minimizing  the  angles  ensures 
that  both  va  and  vb  lie  entirely  in  Q.  It  is  also  easy  to  show  that  the  combination  of  these  two 
segments  resolves  the  reflex  angle  at  v  by  splitting  Q  into  3  polygons  (note  that  in  most  cases,  a  and 
b  will  be  and  v^,  respectively). 

[FIGURE  81 

Figure  8:  Completing  the  convex  decomposition  of  P. 

The  decomposition  algorithm  thus  consists  of  a  recursive  "cutting”  phase  which  relies  on  the  algorithm 
given  for  the  polygon-cutting  theorem.  The  recursion  stops  when  the  polygon  currently  examined  has  fewer 
than  7  notches,  at  which  point  the  procedure  just  described  is  called  upon  to  finish  off  the  decomposition.  We 
observe  that  if  either  the  vertex  v.  or  v.,  say  v.,  in  Theorem  1  is  a  notch,  Le.,  c  =  l,  v.  appears  in  both  of  the 
resulting  polygons  P1  and  P2,  but  is  a  notch  for  at  most  one  of  them.  Therefore  if,  by  extension,  we  let  C(N) 
denote  the  weight  of  P  and  C(AU  (resp.  C(AL)  be  the  weight  of  P.  (resp.  PA  we  can  write: 


(1) 

QiVj)  <  C(W2)  5  2C(A/)/3  +2 

(2) 

Note  that  C(NJ  and  C(N2)  are  actual  weights,  i.e.,  they  count  exactly  the  number  of  notches  in  Px  and  P2, 
as  opposed  to  the  weights  of  P1  and  P2  as  defined  in  Theorem  1,  which  did  not  account  for  the  endpoints  of 
the  splitting  segment  It  is  easy  to  see  that  in  the  worst  case,  we  will  end  up  with  C(N)/6  polygons  with,  each, 
6  notches,  and  the  final  phase  will  use  2  cuts  for  the  resolution  of  each  reflex  angle.  This  will  result  in 
I3C(A0/6  convex  pieces,  which  is  to  be  compared  with  the  minimum  number  of  convex  pieces,  shown  to  be 
always  greater  than  or  equal  to  C(N)/2+ 1  in  [CH80]. 

Next  we  turn  to  the  complexity  of  our  decomposition  algorithm.  While  it  clearly  needs  0(N)  space, 
evaluating  its  run-time  T (N)  calls  for  further  investigation.  If  we  neglect  the  preprocessing  phase  for  the  time 
being,  we  have  the  relations 

T(A0  =  T (JVX)  +  T(N2)  +  O(A0.  ifC(A0  >  6  (3) 

N2  +  N2  =  N  +  2  (4) 


T(A/)  =  O(A0,ifC(A0s6 


(5) 
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Consider  die  recursion  tree,  and  label  each  node  with  the  number  p  of  vertices  in  the  corresponding 
polygon.  At  the  leaves,  we  have  C(P)  <,  6,  while  for  their  ancestors,  C (?)  >  6.  Relations  (3)  and  (5)  show  that 
up  to  within  a  constant  factor,  T(N)  is  equal  to  the  sum  of  the  labels  in  the  tree,  while  from  Relation  (4),  it 
follows  that  if  L(i)  counts  the  sum  of  all  labels  at  level  i,  we  have  L(0)=Afand 

Ui)  s  Ui-l)  + 1 

which  gives  L(i)  +  2i+1.  If  &  is  the  height  of  the  tree,  we  easily  find  that 

L(0)  +  ....  +  L(M)  $Nk  +  2k+\ 

and  since,  from  (2),  we  have  k  =  0(log  C (N))  =  0(log  AO,  including  the  0(  A/log  N )  preprocessing  in  the 
running  time,  we  can  conclude: 

Theorem  7:  In  0(Mog  N)  time  and  with  0(N)  space,  it  is  possible  to  decompose  a  simple 
N- gon  P  into  fewer  than  4.333...XOPT  convex  pieces,  without  introducing  new  vertices,  where 
OPT  is  the  minimum  number  of  convex  pieces  necessary  to  partition  P. 

3.2.  Triangulation 

When  all  the  pieces  of  a  convex  decomposition  are  triangles  and  no  new  vertices  are  introduced,  the 
decomposition  is  called  a  triangulation  of  the  polygon.  An  OfJVIog  N)  algorithm  for  computing  a  triangulation 
of  a  simple  polygon  has  been  given  in  [GJ78],  The  method  requires  the  somewhat  cumbersome  use  of 
A VL -trees  or  similar  logarithmic  search  trees.  This  can  be  avoided  by  using  a  strategy  based  on  the  polygon¬ 
cutting  theorem. 

We  may  choose  to  assign  a  weight  =1  to  each  vertex  of  P  and  apply  the  polygon-cutting  theorem 
recursively,  until  the  polygon  under  consideration  has  fewer  than  7  vertices,  at  which  point  it  is 
straightforward  to  complete  the  triangulation.  We  omit  the  details.  An  alternative  consists  of  computing  a 
convex  decomposition  of  Pas  described  in  the  previous  section,  then  triangulate  each  convex  polygon.  To  do 
so,  pick  any  vertex  of  the  polygon  and  join  it  to  every  other. 

In  both  cases,  a  triangulation  of  P  can  be  explicitly  computed  in  0(Mog  N)  time,  which  matches  the 
performance  of  [GJ78].  We  recall  that  it  is  yet  unknown  whether  kMog  N  is  optimal  for  this  problem. 

Theorem  8:  Using  the  polygon-cutting  theorem,  it  is  possible  to  triangulate  a  simple  iV-gon  in 
0(  A/log  N)  time  and  O(A0  space. 

It  is  shown  in  [CH82]  how  the  additional  information  provided  by  an  arbitrary  triangulation  of  a  simple 
polygon  is  sufficient  to  derive  optimal  algorithms  for  a  number  of  geometric  problems.  Using  the  previous 
result  to  compute  a  triangulation  of  P,  it  is  then  possible  to  determine  the  area  visible  from  any  point  inside  P 
in  linear  time.  Also,  internal  path  problems,  i.c„  problems  involving  the  computation  of  the  shortest  path 
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between  two  points  inside  Pcan  then  be  solved  optimally  [CH82]. 

4.  Conclusions  and  future  research 

The  decomposition  principle  in  geometry  expresses  the  feasibility  of  local  treatments  for  the  solution  of 
general  problems  on  arbitrary  figures.  The  polygon-cutting  theorem  presented  in  this  paper  asserts  the 
applicability  of  this  principle  in  foe  case  of  simple  polygons,  and  by  doing  so,  leads  to  efficient,  simple 
divide-and-conquer  methods  for  solving  a  variety  of  geometric  problems. 

The  merit  of  this  approach  lies  primarily  in  foe  versatility  of  its  applications  as  well  as  in  the  increased 
efficiency  which  it  affords.  The  most  immediate  open  question  is  whether  sorting  foe  vertices  in  preprocessing 
is  indeed  required.  If  not,  foe  algorithm  would  automatically  become  linear.  In  this  paper,  we  have 
deliberately  chosen  simplicity  and  practicality  over  generality  by  restricting  foe  weights  attached  to  the 
vertices  to  take  on  foe  values  0,1.  This  was  motivated  by  foe  fact  that  this  restriction  still  allowed  us  to  apply 
divide-and-conquer  to  an  arbitrarily  chosen  subset  of  vertices,  while  adding  simplicity  to  foe  exposition.  The 
reader  will  observe,  however,  that  it  is  straightforward  to  extend  foe  theorem  to  a  more  general  weight 
function. 

The  applications  mentioned  in  this  work  are  only  a  few  examples  among  a  number  of  other  problems 
which  can  benefit  from  the  polygon-cutting  theorem.  Enlarging  the  list  of  applications  given  here  is  certainly 
a  worthwhile  endeavor. 
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